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Context. Radiative shocks are found in various astrophysical objects and particularly at different stages of stellar evolution. Studying 
radiative shocks, their topology, and thermodynamical properties is therefore a starting point to understanding their physical proper- 
ties. This study has become possible with the development of large laser facilities, which has provided fresh impulse to laboratory 
astrophysics. 

Aims. We present the main characteristics of radiative shocks modeled using cylindrical simulations. We focus our discussion on the 
importance of multi-dimensional radiative-transfer effects on the shock topology and dynamics. 

Methods. We present results obtained with our code HERACLES for conditions corresponding to experiments already performed 
on laser installations. The multi-dimensional hydrodynamic code HERACLES is specially adapted to laboratory astrophysics experi- 
ments and to astrophysical situations where radiation and hydrodynamics are coupled. 

Results. The importance of the ratio of the photon mean free path to the transverse extension of the shock is emphasized. We present 
how it is possible to achieve the stationary limit of these shocks in the laboratory and analyze the angular distribution of the radiative 
flux that may emerge from the walls of the shock tube. 

Conclusions. Implications of these studies for stellar accretion shocks are presented. 

Key words. Hydrodynamics - Radiative transfer - Shock waves - Plasmas - Stars; formation - Methods: numerical 



1. Introduction 

Radiative shocks are shocks in which the structure of the flow is 
affected by radiation, either because the radiative energy density 
cannot be neglected or because the energy transport by radia- 
tion (i.e., the radiative flux) is significant. For densities relevant 
to laboratory astrophysics, the radiative flux becomes important 
at temperature of a few tens of eV. The temperature has to be 
about one order of magnitude higher before the radiative en- 
ergy becomes important. In astrophysics and in particular stel- 
lar physics, the density can be far lower and radiation therefore 
becomes important at lower temperatures. 

A characteristic of radiative shocks is that at high shock 
speeds, the shocked matter is strongly heated and radiates en- 
ergy. Unless the cold upstream gas is completely transparent, it 
absorbs this radiation creating a hot radiative precursor in front 
of the shock where the matter is generally ionized. We focus 
hereafter on shocks that exhibit a developed radiative precursor 
The modeling of these radiative shocks is difficult because 
of the nonlocal coupling between hydrodynamics and radiation, 
of the different scales im portant to radiation and hydrodynamics 
jZerdovich & RaizeiJl9 67: Sincell et al. 1999jj^ and of multi- 
1 (iLevs 



dimensional effects (iLevgnac et al.H2006l: iGonzalez et alj|2006h . 

The topology and evolution of radiative shocks depend on 
several factors, such as the geometry (from ID, as in shock 
tubes, to more complex geometries), the energy deposition (im- 
pulsive, as for radiative blast waves, or continuously driven 
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by a constant piston velocity), the optical depth of the down- 
stream and upstream regions relative to the shock, and the 
possible contribution of a magne tic field. A detailed descrip - 
tion of radiative shocks is given by Zerdovich & Raiz^ (1 19671) . 
iMihalas & Mihaiasl ( [1981 . and lDrakd ilOO^ . 

These strong radiative shocks are found at different stages 
of stellar evolution: from star formation, when the accreting 
gas falls onto the protostellar forming object; duiing the ac- 
cretion - ejection process, which manifests itself in terms of 
the jets of Young Stellar Objects (YSO) and stellar disks, in 
the atmospheres of evolved stars, when shocks dri ve pulsation 
(IChadid et al.1 120001: iFadevev & GilletlfT998ll2000l) : and in the 
last stages of supernovae explosions. 

Spherical or ID symmetry is the usual obvious approxima- 
tion, which prevails in the lack of further information about the 
geometry of these astrophysical shocks. However, if this approx- 
imation is reasonable as a first step for pulsating stars, accretion 
shocks, and supernovae explosions, it may become questionable 
in the case of accretion funnels of classical T-Tauri stars, where 
the matter is supposed to fall from the stellar disk to the stel- 
lar surface by means of channels that follow the magnetic-field 
lines. Accretion in Classical T-Tauri stars is inferred from veil- 
ing in the spectrum and also from signatures in the spectral lines. 
Simulations indicate that the topology of the accretion funnels is 
strongly dependent on the relative orientation of the magnetic- 
field symmetry and rotation axis (Romanova et al. 2003). At the 
base (on the stellar surface), the column may vary from compact 
circular to ring geometries, whereas the upper part (connected to 



2 



Gonzalez, Audit & Stehle: 2D numerical study of radiative shocks 



the disk) may exhibit various fold-like, or more compact topolo- 
gies driven by instabilities. This topic is still the subject of much 
debate and controversy. 

Both the complexity in the structure of the radiative shocks, 
even for a ID approximation, and the lack of angular resolution, 
require accurate simulations coupling radiation transport with 
hydrodynamics. This motivates the study of radiative shocks in 
2D simple cylindrical geometries, to analyze the effect of the 
finite lateral size of the radiative shock, on its structure and lu- 
minosity. One should not use identical simulations to both check 
with and analyze observations. Experimental validation of these 
complex hydrodynamic features is thus required. 

Obviously, the typical dimensions of astrophysical radiative 
shocks ensure that they cannot be studied directly in the labora- 
tory. The study of radiative shocks on Earth thus require different 
gases and physical conditions, which can be recreated by high 
energy installations suc h as laser or pulsed electric installations 
dRemington et al.l2006l) . Radiative sh ock experiments have been 
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performed \yith h i gh power laser s ( Bozier et al. 19861 
Keiter et al.l 12 002'; 'Reighard et al.' '2006: 'Bouquet et al.l 
Gonzalez et al . 2006; Busq uet et al. 2006, 2007). Typically, a 
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200 J laser in about 1 ns can launch a shock at about 60 km/s in 
targets of millimeter size (iBouquet et al.ll2064h filled with xenon 
at pressures of some fractions of bars, whereas supercritical 
shocks at 100 km/s in SiOo aerogel, argon, and xenon have been 
produced at the OMEGA laser (5 kJ). Strong shocks have also 
been generated by a compact pulse-power device CK ondo et al.l 
|2006|) . producing shocks at 45 km/s in xenon with different gas 
pressures of up to 10"^ bars. We see below that, for a fixed shock 
velocity accessible to a laser installation, the radiative regime is 
more easily achieved with a low-density, high-atomic-mass gas, 
which explains why low-density xenon is often used in these 
studies. 

These experiments enable one to examine the physics of 
the radiative precursor (dynamics and density) by visible in- 
terferometry and provide information about the temperature of 
the shock front, by visible emissivity. Studies of the dynam- 
ics o f the precursor of xe non radiative shocks by int erferom- 
etry (iBouquet et al.l |2004|) . and shadowgraphy (Gonza lez et al.l 
|2006|). and of it s topology by instantaneous X-ray imaging 
dVinci et al.l2006l) indicate that multi-dimensional effects can af- 
fect the shock wave, and, in particular, its precursor This was 
attributed to the lateral radiation losses (through the walls of 
the shock tube), which reduce the amount of radiation heatin 
the precursor and thu s affect its structure (iLeygnac et al.ll200 
IGonzalez et al."2006'). 

^eiltyetal. (2000). IShigemori et al.1 (|2000|) 



Edwards et al.l 



2001l) . lCalder et al.l (120021) . and lLaming & Grunl(l2002l) focused 



their studies on the spherical, radiative, blast waves, which are 
similar to those studied in the previous ID cases, in terms of 
the development of radiative precursor, but also exhibit strong 
differences, such as the importance of radiative cooling and the 
development of various instabilities. The study of these experi- 
mental blast waves is linked to the problem of the evolution and 
stability of supernova remnants. 

After the presentation of the typical structure of ID radia- 
tive shocks, we present a numerical study of laboratory radia- 
tive shocks that can be generated by medium-power laser facili- 
ties, such as PALS (Prague, Czech Republic), LULI (Palaiseau, 
France), and similar lasers. The numerical code (HERACLES) 
used in this work is first briefly presented. We then study the ef- 
fect of the duration of the shock driven by the piston (Sect. |4]i. 
In Sect. |5] we consider the effects of multi-dimensional radia- 
tive transfer and their correlation with the width of the shock 



tube. In Sect. |6] we examine the influence of the radiation 
boundary-conditions, which are inferred from the albedo of the 
cell wall, on the topology and propagation of the shock. Finally, 
we present the perspectives and conditions for obtaining an ex- 
perimental stationary radiative shock and present the angular de- 
pendence of the emerging radiative flux of such a shock. 

2. Structure of a 1D radiative shock 

In this section, we present the typical features of a radiative 
shock in terms of ID geometry. When the flow passes through 
the shock, it is heated and it re-emits a fraction of its energy by 
radiation processes. These emitted photons are then absorbed by 
the cold and opaque unshocked gas. This preheated upstream gas 
is called the radiative precursor When the stationary regime is 
reached, a radiative shock is then composed of both a purely hy- 
drodynamical shock (which consists of a discontinuity in the hy- 
drodynamical quantities) and a radiative precursor, whose length 
is directly correlated with both the mean free path of photons in 
the gas and the shock velocity. One usually considers two classes 
of shocks depending on the temperature reached just ahead of the 
hydrodynamical discontinuity. If it is lower than the post-shock 
temperature, the shock is called subcritical, and if it is equal the 
shock is called supercritical. 

By assuming that a strong shock develops in an initially cold 
gas of density po (mass number A) and solving the Rankine- 
Hugoniot relations, one can show that the temperature reached 
in the shock equals 



T = 120.27 K 
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where < Z > is the mean ionization stage in the shocked mate- 



rial, X 



is the ratio of the ionization energy per atom to 



the initial kinetic energy per atom, and ju = the mean molecular 
weight of the non-ionized gas). The photons escaping from the 
shock front are used to heat the cold gas, incoming onto the front, 
to a temperature of T (which is lower th an T). Neglecting the 
gas c ompression in the precursor, one has ( jZel'dovich & Raizen 
Il967l) : 



(2) 



where e - 



^J(y-i? ^ '■^^ internal specific energy (erg/g) after 
taking into account the ionization. 

When T equals T, the shock enters the supercritical regime, 
which is characterized by an extended precursor The corre- 
sponding critical velocity is obtained formally by resolving the 
coupled set of Eqs. ([T]) and (|2|i. This system must be solved nu- 
merically because one needs to take account of the temperature 
dependence of both the internal energ y and ionic charge. U sing 
the atomic -physics model described by Michaut et alj ( |2004|) and 
for xenon of density 10"^^ g/cm^, we found a critical velocity 
around 30 km/s, achieving a shock temperature of 12 eV and 
average ionic charge of 8. In the particular case of a perfect 
gas with 7 = 5/3, and neglecting the ionization-energy term in 
Eqs. ([1]) and one obtains the following approximate value of 



v, > 450 km/s 



< Z > +1 



4/5 



PC) 



1 g/cm-^ 



1/5 



(3) 



which illustrates that the supercritical regime is more easily 
obtained at low density and for high-mass gases. Using ID 
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Fig. 1. Typical profiles of a ID radiative shock: density (upper 
left), velocity (upper right), radiative and gas temperature (lower 
left), and radiative flux (lower right). 



hydro-radiative simulations, lEnsmanl d 19941) found a supercrit- 
ical regime at 20 km/s for hydrogen at 8x10"'" g/cm^ in qualita- 
tive agreement with Eq. |3] 

The precursor extension, L, can be approximated with a ra- 
diative heat-conduction model. One then finds that L is about 
several photon mean-free-paths close to the supercrit ical limit 
jNemchinov & Shuvalovlll980l:lPenzo & Tassartj|1984b . 

A typical ID supercritical shock structure is shown in Fig.[T] 
These academic profiles were obtained using xenon with a con- 
stant ionisation charge of 5, a perfect-gas equatio n-of-state with 
7=1.1, and analytic opacities dBozier et al.lll986l) . The velocity 
of the shock is 65 km/s and the initial density and temperature 
are LSxlO""* g/cnv' and 300 K, respectively. The radiative pre- 
cursor is clearly evident in the radiative flux and temperature 
profiles but, at this moderate shock speed, the density and veloc- 
ity are only slightly affected. 

3. The HERACLES code 

The simulations presented in this paper were completed 
using the three-dime nsional rad iation-hydrodynamics code 
HERACLES ( Gonzalez et al.ll2007h . This code solves the Euler 
equations for hydrodynamics coupled with the moment equa- 
tions of the transfer equation. The original closure relation used 
for the radiative model (Mi model) allows one to be exact 
both in the diffusive and transport limit regimes. This code out- 
put data was verified by analytic models and compared with 
other codes. It can study a wide variety of astrophysical prob- 
lems jGonzalez et al.l |2007|) . Furthermore, it was used to de- 
velop and analyze a l aboratory experiment of radiative shocks 
jGonzalez et alj l2006h . This cross-validation between experi- 
ment and simulation illustrates the relevance of HERACLES to 
both laboratory astrophysics and classical astrophysical situa- 
tions. 

The strength of HERACLES is that it can deal with multi- 
dimensional problems. Until now, radiative shocks have been 
studied by ID geometry models, multi-D models with hydro- 
static fl ows, or by adopting the flux-lim ited diffusion approxi- 
mation ( iBouquet et alJl2064llDrakell2007l and reference therein). 
On the one hand, the multi-D effects (such as lateral losses) can 
determine the structure of the flow, and on the other hand, at the 
foot of the radiative precursor, the reduced flux (ratio of the ra- 



diative flux to the product of the radiative energy and the light 
speed) is nearly one so that we are in the transport limit. In this 
context, HERACLES is of particular relevance to the study of 
multi-D radiative shocks. 

Although all shocks in reality propagate within 3D media, 
we performed 2D axisymmetric simulations. As a first step, we 
considered only the determination of the theoretical influence of 
different parameters on the shock propagation. Furthermore, we 
verified that 2D axisymmetric and 3D Cartesian simulations pro- 
duced identical results (e.g., speeds, positions) if one is careful 
to ensure that the surface-to-volume ratio is equal in both cases . 
Our results are therefore relevant to the case of radiative shocks 
propagating in a rectangular medium, as in most laboratory ex- 
periments. 

Hereafter, all simulations presented are for shocks propa- 
gating in cylindrical cells with the left vertical boundary cor- 
responding to the symmetry axis and the lateral losses located at 
the right vertical boundary. They are filled with xenon, which is 
initially at ambient temperature and a pressure of 0.2 bar (i.e., 
this corresponds to a density of 1.3x10"- ^ g/cm-^). These val- 
ues are typical of lab oratory experiments dBouquet et al J [20041 : 
[Gonzalez et al.ll2006l) where a high molecular weight and a low 
density are chosen to maximize radiative effects. If not speci- 
fied, we use a realistic equation of state and opacities. The equa- 
tion of state of xenon was computed using the OPA-CS code 
of C. Stehle, which uses a screened hydrogenic model and is 
described in iMichaut et al.l ( |2004|) . The opacities were kindly 
provi ded by M. Busquet, who computed them using the STA 
code jBar-Shalom et al.lll989h . The shock is driven at a speed 
of 65 km/s which is a typical of present-day laser experiments. 
The boundary of the cells can have a variable albedo, defined as 
the fraction of the incident radiative flux re-emitted by the wall. 
We demonstrated in previous work dGonzalez et al.l 120061) that 
this parameter has a strong influence on the precursor propaga- 
tion. We therefore performed simulations for albedos between 
0% (fully transparent or fully absorbing boundaries) and 100% 
(total reemission). 

4. Launching effect 

The dynamics of the shock propagation depend strongly on the 
launching procedure. In laser driven shocks, the conversion of 
the laser energy to the piston mechanical energy is performed 
by the shock genera t ed in the pis ton during the laser ablation 
( IBouquet et alj|2004l: lDrakdl2005l) . In this paragraph we com- 
pare the dynamics of the shock for two launching scenarios. In 
the first, the piston moves at a constant velocity. The second 
corresponds to the case where the piston decelerates with time. 
Hereafter, we denote these two driving cases as continuous and 
impulsive respectively. 

Figure |2(a)| shows the difference between a shock with a 
launching phase of the shock that is impulsive and one for which 
the driving is continuous. In both cases, the lateral albedo is 
40%. When the shock driver ("piston") is continuous, the shock 
propagates at constant velocity equal to the initial velocity of 
65 km/s. In the impulsive case, the driver of 65 km/s is continu- 
ous during 0.3 ns and then decreases with a temporal dependence 
that rep roduces qualitatively the velocities of the experimental 
shocks dGonzalez et al.l20"06l) . In that case, the shock decelerates 
as soon as the drive stops because it has little inertia. Regardless 
of the drive duration, one can observe three different stages in the 
precursor dynamics. In the first stage, the precursor moves faster 
than the shock. The photons emitted by the shock propagate in 
the upstream fluid, which is optically thick, where they are ab- 
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(a) Shock and precursor absolute po(-b) Distance between shock and pre- 
sitions cursor 

Fig. 2. Influence of the dynamics of the piston launching on the 
shock dynamics. The shock is driven in xenon at 0.2 bar at room 
temperature. The shock tube is a cylindrical cell with a 0.7 cm 
diameter and albedo of 40%. The full line refers to the impulsive 
lanching whereas the dotted line refers to a shock driven by a 
piston at constant velocity of 65 km/s (cf. Sects. [3]and|4|. 

sorbed. This induces a warm, ionized, radiative precursor in the 
upstream part of the shock. In a second phase, the precursor de- 
celerates (it can be even slower than the shock front). Finally, it 
reaches a constant speed equal to that of the shock front. This is 
the stationary regime. 

These three stages are clearly seen in Fig. |2(b)| which shows 
the distance separating the shock from the precursor leading- 
edge versus time. Since the precursor goes faster than the shock, 
this distance increases to a maximum value. Then, as the precur- 
sor slows down, this distance decreases until reaching an asymp- 
totic value (stationary limit), when both the shock and precursor 
leading-edge travel at the same speed. 

In a laboratory experiment, the duration of the driving is not 
a parameter that one can vary simply. It is easier to modify the 
shock dynamics by changing the radiative lateral losses, which 
can be achieved by varying the width of the canal or by mod- 
ifying the albedo of the cells walls. This albedo can be spe- 
cially adapted by choosing properly the material of the cell or 
by choosing special coating on these walls. 

In the following sections, we examine the influence of these 
parameters on the dynamics of the shock precursor and in par- 
ticular on the delay before the stationary regime. 

5. Effect of a 2D radiation field 

We define R to be the ratio of the transverse width 1 to the pho- 
ton mean-free-path A in the unshocked medium (A ^ I mm). In 
this section, we present the influence of R on the radiative shock 
structure and evolution. In our bi-dimensional and axisymmet- 
ric simulations, the gas is driven at a speed of 35 km/s, constant 
both in time and across the cell. Since we impose lateral reflect- 
ing boundaries on the hydrodynamics, the shock would remain 
planar without accounting for the radiative transfer. It would be 
transversally homogeneous, and temperature and density would 
vary only with y-position (along the canal axis). However, when 
accounting for radiative transfer with lateral losses through par- 
tially reflective walls, 2D effects can be observed. Density and 
temperature vary with x-position (perpendicular to the axis) and, 
in some conditions, curvature of the shock front may even be ob- 
served. The influence of the radiation on the shock structure is 
directly correlated with the possibility of photon escape from the 
medium, which can be estimated from the dimensionless num- 
ber R=l/A. 

The importance of radiative losses on the shock structure, 
which depends strongly on R, has not been considered in de- 



tail by many laboratory experiments nor in the study of radia- 
tive shocks encountered in astrophysical objects. For example, 
the mean free path in xenon at 300 K and lO^^* g/cm^, which 
is typical of laboratory experiments, is about 0. 1 mm whereas 
the canal width is about 1 mm. For a typical stellar atmosphere 
(T=5000 K and p=5xlO"" g/cnr'), the typical mean free path is 
4x10'^ cm and decreases rapidly when the temperature increases 
dSeaton et al. 1994) whereas the radius of a classical T Tauri star 
is about 2x10" cm. 

Figure [3] shows bi-dimensional maps at the same time for 
three different canal widths (i.e., three values of R). The shock 
position is marked by the jump in density and by the temperature 
peak. The precursor is located in the lower part of the figure 
below the shock. Its extension is particularly visible in the maps 
of temperature and radiative flux. 

If R is small (cf. Fig. |3(a)[ ), the photons don't interact with 
the gas and they escape freely from the system. As the radia- 
tion escapes laterally, the shock loses some part of its energy, 
it decelerates rapidly, and the precursor remains thin. Since the 
photons emitted at any point of the downstream medium es- 
cape freely, the shock remains planar. In contrast, if R is large 
(cf. Fig. |3(c)[ ), the photons cannot escape because they are re- 
absorbed over a short distance. The shock travels faster than 
in the first case. In the limit where R tends to infinity, we re- 
cover the uni-dimensional case with a planar shock except for a 
small layer along the walls, similar to a boundary layer. In the 
intermediate case (cf. Fig. |3(b)| i, the shock is curved, since only 
photons close to the walls escape, and photons emitted on the 
axis are trapped. The shock then decelerates more close to the 
walls than on the axis, which tends to bend the shock. Shock cur- 
vatures hav e already been observed experimentall y by imaging 
techniques jVinci et alj2006l:lReighard et alj2006l) . although the 
dynamical origin of the curvature can be of many different na- 
tures. 

6. Effect of radiative losses 

In the previous section, we showed that the ratio R of the pho- 
ton mean free path to the canal width was a key parameter for 
the radiative losses. In this section, we study in more detail the 
influence of the lateral radiative losses on the shock and precur- 
sor dynamics. We therefore consider a canal of fixed width but 
with walls of variable albedo a, which varies between (no wall 
reflexion) and 1 (purely reemitting walls). 

This wall albedo is a crucial point in the laboratory experi- 
ments. Varying the element used for the cells or covering the cell 
walls with layers of different elements (e.g., gold, aluminium), 
one can test the radiative-loss effects. 

We study the influence of the albedo on the dynamical prop- 
erties of the shock at short and long times. To illustrate the ef- 
fect, we have chosen to focus the analysis on the case of im- 
pulsive launching. We concentrate our study on the positions of 
the shock front and of the leading edge of the precursor (defined 
as the point, on the cell symmetry axis, where the temperature 
reaches the threshold value of 5 eV). The value of the temper- 
ature threshold influences the extension of the precursor only 
slightly. 

Figure ID shows the shock dynamics obtained with an albedo 
of 0%, 10%, 20%, 30%, 40%, 50%, 60%, 70%, 80%, 90%, and 
100%. One can see that the wall albedo has almost no effect on 
the shock velocity, at least for the evolution times considered 
here. However, the precursor is strongly affected by the albedo. 
As can be seen in Fig. |4(b)| for small values of the albedo (for 
important lateral losses), when most of the photons escape, the 
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Fig. 3. At a same instant, snapshots of the maps of density, tem- 
perature, radial, and axial radiative flux showing the influence of 
the canal width compared to the photon mean free path (R=0.5 
top panel, R=l middle panel and R=10 bottom panel). Case 
of xenon initially at 1.3x10^^ g/cm^ and 300 K, using perfect- 
gas equation-of-state {y-5/3 and mean molecular weight fi-20 
equivalent to a mea n ionization stage of 5), and analytic opacities 
jBozier et al 1 119861) . The shock propagates from top to bottom. 
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Fig. 4. (a) : Influence of the walls albedo upon the shock (dashed 
line) and precursor (solid line) positions, (b) Influence of the 
albedo upon the distance between the precursor and the shock. 
Comparison for an albedo of 0%, 40%, 50%, 60%, 70%, 80%, 
90%, and 100% (from bottom to top). 



precursor extension remains small and its deceleration increases. 
This dependence is non-linear (cf. Sect. |2]below for the depen- 
dence of the time needed to reach the stationary limit on the 
albedo). The extension of the precursor relative to the shock 
front reduces slowly with time, as a consequence of the weak- 
ening of the shock and the slow deceleration of the piston. This 
evolution at later times differs from the case where the piston 
moves at a constant velocity, as can be seen from Fig. [21 which 
converges at long times towards a stationary limit, where precur- 
sor and shock front move at the same velocity. At early times, the 
precursor front moves far more rapidly than the shock front, and 
we expect that this regime is close to the case of a Marshak wave 
(a radiation-driven thermal wave without coupling to hydrody- 
namics). For t below 5 ns, and a=l, one thus recovers the well- 
known ID-Marshak-wave time-dependence of yfi. For smaller 
values of the albedo, the structure of the 2D-bent Marshak waves 
is more c omplex, even supposing a con stant opacity, as one can 
see from [Hurricane & Hammerl ( |2006|) . who studied the bend- 
ing and slowing of the front by radiative losses. However, even 
if similarities exist between the dynamics and topologies of a 
Marshak wave and those of a radiative precursor at early times. 
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the two processes differ physically, especially at late times where 
the stationary limit is reached in one process and not in the other. 

7. Toward stationary shock 

In the case of a shock with continuous driving, the typical evo- 
lution of the radiative shocks exhibits, as in the previous case, a 
first phase of development of the precursor, where the ionization- 
front velocity exceeds the shock velocity. After this phase, the 
ionization-front velocity decelerates until it becomes equal to 
the front velocity. At this time, the stationary limit is reached. 
The corresponding time at which this stationary radiative shock 
is formed, f 5,3,, is of the order of the ratio of the precursor exten- 
sion to the shock velocity. 

The stationary shock limit is usually assumed for astrophys- 
ical flows, because f,,,,, is less than the typical hydrodynamical 
time, t/jyij- For instance, in the case of pulsating stars, f/,,.^/, which 
is roughly equal to the ratio of the atmosphere width to the shock 
velocity, varies between 1500 s for RR Lyrae and 10^ s for cooler 
giants such as RV Tauri, whereas t stm is of the order of 50 s for 
a shock velocity around 10 km/s jFokin et al.ll2004l) . In accre- 
tion shocks along magnetospheric columns, this time is approxi- 
mately equal to the ratio of the column length to the shock veloc- 
ity, or 10"*- 10'' s for classical T Tauri stars (of about 3 solar radius 
and 300 km/s). Thus, in the absence of any shorter timescale of 
hydrodynamical motion (e.g., instabilities), these astrophysical 
shock waves can be supposed to be stationary from the radiative 
point of view. 

Associated with the ID approximation, this stationary limit 
also offers the advantage of simplicity in the observational di- 
agnostics w ith the help of publishe d stationary-shock structure 
tabulations (iFadevev & Gilletll2000l) . 

In the experimental cases, the piston velocity decreases 
slowly with time. However, as in the academic case of contin- 
uously sustained shock, the precursor extension relative to the 
shock front increases at early times until reaching a maximum 
value. We have chosen to define the time interval to form a sta- 
tionary shock, tsiai, as the time of maximum extension of the pre- 
cursor because it unambiguously delimits the short accelerating 
and long decelerating phases of the precursor dynamics. 

Using this definition, and in the launching case, whic h is 
close to that of previous experiments dGonzalez et al.ll2006l) . we 
find tstat = 5.6 ns, 6.1 ns, 7.1 ns, 7.6 ns, 8.6 ns, 10.6 ns, 13.1 ns, 
18.7 ns, 32.3 ns, and 409 ns for wall albedos of 10, 20, 30, 40, 
50, 60, 70, 80, 90, and 100 percent, respectively (cf. Fig. |5]l. 
Thus, high lateral radiation losses (small albedo) not only decel- 
erate the radiative precursor, but also strongly decrease the time 
to reach the stationary limit. From these numerical simulations, 
we derive an analytical fit 



1000 



(fi - fo)(l -ay + to 



(4) 



where a is the wall albedo, to and fi are the stationary times for 
a=0 and a=l, respectively, and /? is a constant. Their values are 
to = 5.6 ns, ti = 409 ns, and = 0.84. 

In all of these simulations, the shock position is insensitive 
to the albedo of the tube walls. However, even if the shock speed 
is constant, the density and temperature of the post-shock mat- 
ter are sensitive to the radiative losses. Higher losses (smaller 
albedo) result in a higher density and a lower temperature. As 
a consequence, the precursor extension tends to be lower in the 
case of higher losses because the emissivity of the shock front is 
proportional to T"* and also because a large fraction of the pho- 
tons escapes without contributing to the precursor extension. 




40 60 
Albedo (%) 



100 



Fig. 5. Influence of the albedo upon the time interval to form 
a stationary radiative shock, tstat (points: numerical simulations 
results, line: fit). 



Until now, we have emphasized the strong dependence of 
the stationary-shock structure on its radiative parameters (e.g., 
ratio R of the transverse width to the photon mean free path, 
and the percentage of radiative losses at the walls of the shock 
tube), or equivalently to departures from the ID approximation. 
Another interesting quantity is the radiation flux, since it is used, 
for instance, as a tracer of accretion rate. In Fig.|6] we plot, for 
different tube-wall albedos, the variations in the luminosity as 
a function of the angle of observation 9 [0 = corresponds to 
a forward flux, 9 - n/l to a perpendicular flux, and 9 - n to 
a backward flux). The plotted profiles correspond to the lumi- 
nosity integrated from the beginning of the radiative precursor 
to the shock front at time t=50ns. The integral is evaluated by 
calculating the photon angular distribution function, which is 
known from the Ml model, at each point along the cell wall. 
We plot a normalized luminosity because we are interested only 
in the angular distribution and not in the absolute value, which is 
more difficult to compare. One can see that with increasing lat- 
eral losses the angle corresponding to the maximum of the lumi- 
nosity increases and the distribution sharpens around this angle. 
This is due to the fact that when the losses are high, the radia- 
tive reduced flux is higher, and the luminosity tends to become 
almost mono-directional (in the limiting case of a unit-reduced 
flux) rather than an isotropic Planck function (in the limiting 
case of diffusion). The luminosity is much more directional and 
the dispersion around the preferred direction is lower. Figure Q 
shows that this angle is approximately a linear function of the 
lateral albedo of the tube walls. All these points illustrate the 
strong anisotropy of the flux in the radiative shock. Although il- 
lustrated for a given time, the conclusions about the influence of 
the albedo on the peak of the angular distribution and on the dis- 
persion about this value are of a quite general character during 
the shock propagation. 



8. Conclusion 

We have illustrated how radiation influences the radiative-shock 
topology and structure. These effects can be observed with ex- 
perimental results and used to interpret observational data. 

We have in particular shown that lateral radiative losses in- 
duce a curvature in the shock front and a strong shortening of 
the radiative precursor. For laboratory experiments, it means that 
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Fig. 6. Normalized luminosities for different values of albedo 
(80% in yellow solid line, 70% in red dotted line, 50% in green 
dashed line, 20% in blue dashed-dotted line and 0% in black 
dashed-double-dotted line). 
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Fig. 7. Angle corresponding to the maximum of the luminosity 
as a function of the albedo (points: numerical simulation results, 
Une: fit). 



the tube geometry and the opacities play a key role in the shock 
propagation and topology. In astrophysical cases, this suggests 
that the link between accretion rate and photometric observa- 
tion is strongly dependent on the gas conditions and the possible 
presence of dust. 

Lateral radiative losses through the walls of the tube lead 
to a weakening of the shock, and a reduction in the precursor 
extension and temperature. They also strongly reduce the time 
needed to reach the stationary limit. 

The radiative flux that emerges from the tube boundaries 
evolves from an isotropic situation in the case of low radiation 
losses to a more anisotropic distribution in the case of strong 
radiative losses. This effect may have a strong impact on photo- 
metric observations of anisotropic accretion flows, such as in the 
case of T Tauri star accretion through funnels. 

High-energy, high-power installations, such as lasers, allow 
generation of strong radiative shocks with extended precursors. 
Most studies have focused on the global characterization of the 
radiative precursor. Only a few have attempted to study the shock 
front, which requires X-ray diagnostics of high temporal and 
spatial resolution and is thus very challenging. An important re- 



sult of our study is the reduction with radiative losses in the time 
needed to reach the stationary limit. Although some experiments 
have been close to reaching this limit (Gonzalez et al. 2006), this 
has not been studied until now experimentally. 

Shock velocities of 60 - 80 km/s are accessible with kJ 
class lasers. Usin g more powerful las er-energy installations 
such as the NIF ('Haynam et al.' '2007^ at Livermore (USA) 
or the LMJ dBettin g er & Decroisett e 1999; Besnard 2Q(M 
lEbrardt & Chapul 20081) at Bordeaux (France) will allow us to 
explore the effects of the radiation pressure on the structure of 
the shock wave. Using xenon at densities of the order of 5x10""* 
g/cm^, such effects start to be visible fo r shock velocities of the 
order of 200 km/s (iMichaut et al.ll2004l) . 

Another interesting topic that could be addressed is the even- 
tual development of radiative instabilities. These "thermal" in- 
stabilities have been studied theoretically and numerically for 
optically thin radiative shocks, in the context of, for exam- 
ple, accretio n onto white dwarfs, and coll i ding radiative flows 
(iFaiia [1971 IChevalier & Imamural [T9p iBertschingeil [19861: 



Rvu & Vi shniac"l987': I Walder & Folinill 19961: iLaming & GrunI 
2002|; [Mignone 2005). Such an instability was observed exper- 
imentally in 3D blast- wave devel opment over a few hundreds 
of nanoseconds jGrun et al.lll99ll) . although oth er experiments 
failed to reproduce it on a shorter timescale jEdwards et alj 
1200 Ih . The development of these instabilities in an optically- 
thicker radiative shock (with a radiative precursor) has, to our 
knowledge, not been studied theoretically. The multidimensional 
effects, by deforming the front and allowing a transition towards 
an optically thin case, could be favorable to the development of 
such instabilities, and this will be the subject of future studies. 

These laboratory studies allow us to test the adequacy of 
modern codes in dealing with astrophysical, hypersonic radia- 
tive flows. The HERACLES 3D radiation hydrodynamics code 
is well adapted to the study of these anisotropic radiative flows. 
It allows us to test together the cases of both optically thick 
and optically thin regimes. It has been already used in model- 
ing the propagation o f astrophysical ra diative jets in the ISM 
dGonzalez et al.l 120061 : iGonzalej |2006[) . Future improvements 
will include the inclusions of multigroup radiative-transport 
(which may influence the fine structure of the shock precursor), 
the decoupling of ionic and electronic temperatures (which af- 
fect the shock front), and non-LTE effects, which play a role in 
the optically thin regions of the shock structure. 
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